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Abstract 

A new theory on gas-liquid phase transition is given. The new idea is that the total intermolec- 
ular potential energy for a classical system in equilibrium is relative with the average distance of 
molecules. A new space homogeneity assumption is postulated, a new formulation — the mean 
distance expansion, instead of Mayer's Cluster Expansion, is introduced, a new explanation on 
liquid-gas phase transition is given, and the physical quantities in a system of coexistent vapor and 
liquid — the densities of the coexistent vapor and liquid, the vapor pressure and the latent heat, 
are calculated. 
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I. INTRODUCTION 



In my recent paper a new theory to discuss first order gas-liquid phase transition in 
statistical mechanics is given. The basic idea is the following: 1) A new Space Homogene- 
ity Assumption (SHA) is postulated: In a classical system in equilibrium with volume V, 
particle number A^, when gravity and boundary effects are omitted, each particle occupies 
volume V/N averagely. The SHA is quite different from Eyring's assumption(it is referred 
as EHSA SHA later) j^: the number of the particles in a volume d^r is N/VAirr'^dr. We 
will demonstrate in detail later that the SHA, in contrast with EHSA SHA, will exhibit the 
correlation between the mean distance of molecules and the partition function. 2) In order to 
calculate the intermolecular potential and the partition function for a classical system, a new 
formulation — mean distance expansion, instead of Mayer's cluster expansion j^, is intro- 
duced, which is characterized by the intermolecular mean distance r 

as a parameter, appearing in the expressions of the intermolecular potentials and the parti- 
tion functions. This implies that the intermolecular potentials and the partition functions 
are dependent on the intermolecular mean distance. Thus the derivatives of partition func- 
tions (some thermodynamic functions) will be relative to the intermolecular mean distance. 
3) A new explanation on first order gas-liquid phase transition in statistical mechanics is 
given: First order gas-liquid phase transition does not correspond to singularities of ther- 
modynamic functions. For a given temperature T < Tc and some pressure P, there are 
two densities of particle number Ug and rii (densities of vapor and liquid) that satisfy the 
equations P{T,ng) = P{T,ni) and fi{T,ng) = fi{T,ni). In each branch thermodynamic 
functions, such as fig (T,ng) or /i; {T,ni), are analytic. 4) The new theory tells us how to 
calculate the physical quantities in a system of coexistent vapor and liquid: densities of 
coexistent vapor and liquid, vapor pressure and latent heat, etc. 

In this paper, we will give some concrete models, as examples, to exhibit how the new 
theory works. In section two, starting from the SHA and the mean distance expansion, and 
by using Sutherland potential, we obtain partition function of VDW fluid and coexistence 
of two phases is discussed. Densities of vapor and liquid, vapor pressure and latent heat for 
VDW fluid are also calculated. In section three, modified VDW fluid is discussed and the 
better results are obtained. In section four, the differences between new theory and Mayer's 
cluster expansion are identified. 
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II. MEAN DISTANCE EXPANSION OF THE VDW FLUID 



For a classical system in equilibrium, we consider the pair interaction which is only 
relative to the distance of particles. The potential energy between particle i and particle 
j is u (fij), where is the distance of these two particles. We assume that the potential 
energy u {vij) is shared by particle i and particle j together and is divided to these two 
particles equally. According to the SHA, particle i occupies a volume V/N = 1/n, where 
n is the particle number density. For convenience, suppose that this volume is a sphere 
with radius ri = (3/47rn)^''^, and particle i is located at its center. We then draw many 
spherical surfaces centered by particle i outside the sphere, so that the volume of every 
shell formed by neighboring two spherical surfaces is V/N, and we will have one particle 
per shell averagely. The radius rj of the j'th sphere is determined by 47rr|/3 = j/n. Then 
intermolecular potential energy of particle i can be expressed as 

N-l 

and the total potential energy of the system is U = NUi = ^ u {rj). Thus the partition 

i 

function is given by 

/exp [-^U) dq = ]^exp jdl (2) 



3N 



where A is thermal wavelength. If the intermolecular potential is the Sutherland potential 

u{r) = { (3) 
-e{a/ry, r > a. 

We obtain the partition function of VDW fluid 

with a = Yli ^^lf~' ^ ~ 1^^"^' ^^"^ ^^^^ volume f dq = (V — Nb) . 



A. The Equation of State 



The VDW equation of state is obtained as follows 



P = knT 



an (5) 
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As well known, we have the critical data 

1 

'^c = — r 
36 



The reduced equation of state is 



Pc 



rirknTr, 8 



"i — n 

with reduced quantities P* = P/Pc,T* = T/Tc and n* = n/n, 



(6) 



ksT, = ^ (7) 
276 ^ ' 



(9) 



P* = - 3n*2 (10) 



B. Chemical Potential 



Chemical potential is expressed as follows 
d\nZ' 

- V,T 



-ksT 



dN 



n 



1 — nb 



-2an + ksT In 
We find that the solutions of equations (12) and (13) 

0, 



nksTb SksT 2TimkBT 
+ ; z — In ■ 



nh 



dn 



0. 



are the same as the formulae (6)-(8). Letting 

7^ 



/i 3 IQirma 



and substituting Eq. (11) into (14), we have 



■Jl{T*,n* 



In 



n 



n 



-InT*-, 

3-n* 2 3-n* AT*' 

which is a function of reduced quantities only. Now we can define the reduced chemical 
potential Ji* = Ji/jic, and draw its isotherms as shown in Fig.l. 



9n* 



(12) 



(13) 



(14) 



(15) 
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FIG. 1: the reduced chemical potential at different temperature for VDW fluid 



C. Densities of Coexistent Vapor and Liquid 



For a coexistent gas-liquid system, we have the conditions for two phases to exist in 
equihbrium 



P* {T\n*g) = P* (T*, 



,ni), 



(16) 



r {T\n;)=r {T 
From equations (10) and (16), we have 

8T* 



(3-n;) {3-n*) 
From the equations (14) and (17), we have 



(17) 



(18) 



4X* 



n 



9 [nl - n:) . 



(19) 



(3-n;)(3-r2;)_ 

The solutions of Eq. (18) and Eq. (19) given by numerical method are plotted in Fig.2. 
This curve is equivalent to that given by Maxwell rule in thermodynamics. 

The curve in Fig.2 can be described by the empirical formulae (firstly introduced by Wei 
Wang ) 



71*1 = 1 + 0.497 (1 - T 



*\ 1.086 



0.56(1 -T* 



,1.53 



2fl -T' 



xO.5 



(20) 



n*g = l + 0.497 (1 - T*f-^^^ + 0.56 (1 - T*f-^^ - 2 (1 - T*f-^ 



(21) 
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FIG. 2: Densities of coexisting vapor and liquid: blue — the experimental curve from Guggenheim 
green — VDW fluid, red — MVDW fluid. Both of red and green are the theoretical results in this paper. 



or 



^ = -0.56 (1 - T*f-^^ + 2 (1 - T*f^ , (22) 



' ^ ' = 1 + 0.497(1 -T*)'■°^^ (23) 

Formula (22) suggests the critical exponent (3 = 0.5 when T* ^ 1. Formula (23) states 
the well-known law of the rectilinear diameter. Comparing with the numerical values of the 
solutions of Eq.(18) and Eq.(19), the inaccuracy of the formula (20) is less than 0.1% when 
T* > 0.6, and the inaccuracy of the formula (21) is less than 0.1% when T* > 0.71. The 
percentage of inaccuracy increases as the temperature drops below T* = 0.71. 

D. Vapor Pressure 

The reduced vapor pressure is expressed as 

- 1^ - '"f - - (^^' 

Substituting Eq. (19) into Eq. (24), we obtain another expression of vapor pressure 

p. ^ 3nf exp {nt - n;) - ^ {nf - nf)] - nf 
l-exp[^{n*-n;)-^{nf-nf)] 
Fig. 3 illustrates the relation between vapor pressure and its corresponding temperature. 
And this relation can be expressed by the empirical formula 

In P* = 3.5204 - 3.5660/r*. (26) 
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FIG. 3: Comparison between the experimental and the theoretical results of vapor pressure. Blue- the 
experimental curve from Guggenheim [gj , green — VDW fluid, red — MVDW fluid. 



E. Latent Heat 

The free energy of the system can be written as 

F = keT In N\ - -NkeT In aNn - NkeT In {V - Nb) . (27) 

The internal energy reads 

= ^NkeT - aNn. (28) 

and then we obtain the difference between the internal energies per particle in two phases: 

At/, = ^ - ^ = a (n, - rig) = (r,) - (r,) , (29) 

where Ui (fg) and Ui (n) are the potential energies per particle in vapor and liquid deter- 
mined by equations (1), (3), (20) and (21) correspondingly, which are only dependent on 
mean distances in the two phases. Multiplying the Avogadro's constant N^, we obtain the 
difference between the molar internal energies in the two phases 

AU = NAa{ni-ng). (30) 
The first law of thermodynamics is expressed as 

AU = AQ- P{vg - vi) =L- P{vg - vi), (31) 



E 



d_ 
df 



F 
T 
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FIG. 4: Comparison between the experimental and the theoretical results of latent heat. Blue- the ex- 



perimental curve of Ar from green — VDW fluid, red — MVDW fluid from Eq.(48) with M = 2.83 and 
S = -0.49 



where Vg and vi are the molar volume of vapor and liquid respectively. And we obtain the 
latent heat for VDW fluid 



NAa{ni - Ug) + P{vg - vi) = Na 



1 1 

a{ni - Ug) + P{ 



n„ 



ni 



(32) 



which is equivalent to Clapeyron's equation(see Appendix A). We define L = L/PcVc and L 
can be written as 



L = 3(n;-<)+P*( 



1 



1 



(33) 



which is a function of reduced quantities p*(T*) and n*{T*) only. Fig.4 shows the relation 
between latent heat L and temperature T*. The blue curve is drawn from experimental data 
of argon The green one is from Eq. (33). 



III. MODIFICATION TO THE VAN DER WAALS EQUATION OF STATE 

In 1945, E.A.Guggenheim collected experimental data on densities of coexisting vapor- 
liquid for eight substances and plotted on a curve. He found that the curve could be 
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represented by the empirical formulae 



^ = 1 + 0.75(1 - ^) + 1.75(1 - ^y/'^ (34) 

^ = 1 + 0.75(1 -|)- 1.75(1 -|)V3 (35) 
It means that the critical exponent 3=1/3 when T* ^ 1 for real erases. 

n 

It was also E.A.Guggenheim |6| who provided the empirical formula on the equilibrium 
vapor pressure 

5 SI 

lnP* = 5.29-^ (36) 

Comparing with Guggenheim's experimental data (34)-(36), we find that the formulae 
(20), (21) and (26), like in the case of the latent heat (33), are right qualitatively, but not 
right quantitatively. The reason is that the Sutherland potential Eq. (3) corresponding to 
the VDW fluid does not reflect the actual interaction. In order to obtain the better results 
and to keep a form as simple as VDW equation, we introduce the effective potential 

+00, r < cr, 
U{r) = { 37 
^ -£(a/r)^^T^, r>a. 

where M and S are constants. In the Mayer's cluster expansion, to keep convergence we 
have to let M ^ 3. But in the mean distance expansion there is no such restriction. 
Following the standard procedure like in section two, we have the partition function 

^ exp{paT^Nn^){V - Nbf (38) 



N\X 



3N 



Eq. (10) reads 



p* 



12(3 + M)n*T*(i-^) 6 + M , a/, 



M(6 + M - Mn*) M 



-n 



(39) 



with 



M , , 

= (6TM)fe 
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(41) 



Pr = — -?^-—^Tf 



when M = 3, S = 0, the 
the Berthlot's equation 

Eq. (11) and Eq. (15) read 



q. (39) is VDW equation, when M = 3, 5 = -1, the Eq. (39) is 



1 + y 1 an^T^ 



kBTln 



n 



nb 



+ ^ , ^In TTT^- (43) 



1 — nb 



/i = In 



n*M 



n*M 



6 + M-n*M 6 + M-n*M 2 



-InT* 



n 



=fT*(5-i)(6 + M)^ 



12M 



Eq. (24) reads 

rp*{l-S) 

Eq. (25) reads 

rp*{l-S) 

Eqs. (32) and (33) read 



(6 + A/ - Mn*g) (6 + Af - Mnl 



12(3 + M) (n! -n;) 



{6 + Mf{ng~-n*~] 
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nj(6+M-Mn;*) M{e+M){n*-n*) 

n* (6+M-Mnj) (6+M-Mn*) (e+M-Mri;*) 



M 



M M 



L = N^iaT^il - 5)(V - < ) + P(- - 1)] 

Ug ni 



(44) 



(45) 



(46) 



(47) 



' +P*( 



1 



1 



(48) 



The solutions of Eqs. (44) and (45) with Af = 2.83 and S = -0.49 are plotted in Fig.2, the 
corresponding vapor pressure and the latent heat are shown in Fig.3 and Fig.4. All of them 
are better than the results of VDW fluid. 

As shown in Fig.2, the curve of coexisting vapor-liquid densities for the modified VDW 
fluid locates between the VDW curve and Guggenheim's experimental curve, indicating that 
the critical exponent (3 for the modified VDW fluid values between 1/2 and 1/3. How can 
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FIG. 5: The gas-liquid coexistent curve(red) obtained by MVDW fluid with S(T) and M(T) given by 
formulae (48) and (49) almost exactly coincide with Guggenheim's experimental curve(blue), with (3 = 1/3. 



we get /3 = 1/3 so that it coincides with experimental data? By letting M and S in Eqs.(44) 
and (45) as the functions of T* as follows 



M{T* 



(T* - 0.69)2 



0.30862 



0.998 -T* 



IT* - 0.9981 



+ 3.865 



(49) 



5(T* 



-(4 - AT 



*\-0.45 



0.92)(T* -0.7) 



0.38 



(50) 



|T* - 0.7| 

we obtain densities of vapor and liquid and n^*, which almost exactly coincide with 
Guggenheim's experimental curve, with /? = 1/3, as shown in Fig.5. 



IV. DISCUSSION 



The partition function of VDW fluid was derived by J.E.Mayer and M.G.Mayer in 1940s 
from their cluster expansion, and by H.Eyring, D.Henderson, etc. in 1964 from their EHSA 
SHA. What are the differences between this paper and their works? We will find out in 
detail in this section. 



A. The New Idea in this Work 

The new idea in this work is that the total intermolecular potential energy for a classical 
system in equilibrium is relative to the average distance r of molecules. It is also the key 
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difference between this work and Mayer's cluster expansion theory. To demonstrate this, we 
draw an isotherm of the reduced pressure p* and the reduced potential energies per particle 
U* = Ui{n)/Ui{nc) against the reduced densities n*, as shown in Fig. 6. When the system 
runs from a to g, or from / to d, the reduced potential energies per particle will vary as n* 
changes, which is completely different from Mayer's work. 

B. Existence of Two Phases 

In a homogeneous system with a temperature T < and the corresponding vapor 
pressure P, when the number of particle is unknown, the phase in which this system exists, 
for example liquid phase or vapor phase, is uncertain. These different fluid phases correspond 
to different solutions of Eqs.(16 - 17), which form two branches of solutions of Eqs.(16 - 
17). 

C. Singularities on Isotherm 

In Fig. 6 on the isotherm a — g — I — d there exist singularities at the point g or /. The 
function P = P{T*, n*) is not different iable at (7 or / on the isotherm a — g — I — d. This is 
due to the fact that the segments a — g and g — l in Fig.6 are described by different partition 
functions respectively. The segment a — or / — d is on an isotherm of a homogeneous 
system, which is described by partition function (4). The segment (7 — / in Fig.6 describes 
an isotherm for an equilibrium coexisting vapor-liquid system that is not a homogeneous 
system and cannot be described by partition function (4). In Appendix B we give a partition 
function for an equilibrium coexisting vapor-liquid system, which corresponds to the segment 
9-1. 

D. Nature of Latent Heat 

The calculation about latent heat is one of the most important parts in this paper, 
which provides a powerful support to the SHA — the basic assumption in this work, and 
distinguishes this work from previous works, for instance, Mayer's cluster expansion and 
the EHSA SHA, evidently. The average distance of gas molecules is larger than that of 
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FIG. 6: The relation between the potential energy per particle U* and densities n*.The curve a-g-b-e-l-d 
is an isotherm of pressure P*. When pressure runs from a to g, the densities vary from n* to n*, the 
corresponding potential energy per particle U* change from U*(n*) to U*(n*).The gap U*i — U*g between 
two phase is shown 

liquid molecules. The interaction of gas molecules and the interaction of liquid molecules 
can be described by the same interaction potential, such as Eq. (2). The difference between 
the mean distances of molecules in the two phases results in the difference of molecular 
potential energies in the two phases. When a liquid molecule vapors and becomes a gas 
molecule, there is an interaction potential energy gap between the molecules in the two 
phases Af/j = Ui(f^) — Ui{rj) = a{ni — Ug). On the other hand, the work done by the change 
between the volumes per mole in the two phases is P{vg — vi). The sum of the interaction 
potential energy gap and the work is equal to latent heat. This is the true physical meaning 
of latent heat — what Mayer's cluster expansion and EHSA SHA have failed to point out. 

E. Analytic Properties of Thermodynamic Function 

From the discussion about the liquid-gas transition in this paper, we know that the 
gases and the liquids appear in different branches of diagram P — n and diagram fi — n, 
but yet surface tension is not taken into account. When surface tension is considered, 
supersaturation and supercooling may occur. The isotherm will follow a — g — hoid — l — e 
in Fig. 6. On the segments g — b or I — e, the system is not in a stable equilibrium, but still is 
homogeneous, and can be described by partition function (3). In these cases, thermodynamic 
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functions are analytic in each branch, even at the first order transition points g or I. 

When T T^^ we have Ug^ric and ni^Uc, the critical phenomena will appear. The 
coexistence of two phases given in this paper will play a very important role in critical 
phenomenon. We will discuss it in another paper. 
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APPENDIX A: EQUIVALENCE BETWEEN EQ. (31) AND THE CLAPEYRON 
EQUATION 

Eq. (31) can be written as 

L = Ug-Ui + Piv, - vi) (Al) 
Applying thermodynamic relation 

U = fi + Ts - Pv (A2) 

and phase conditions 

P, = Pi 



fig = (A3) 

we have 

Ug-Ul= T{Sg - Si) - P{Vg - V l) (A4) 

then Eq. (Al) is 

L = T{sg - Si) = T^^^ivg - vi) (A5) 

According to Maxwell relations and the characters of two phases in equilibrium j8| 

,ds ^ , dP , , . , 

14 



.ds^ Sg- Sl 



It 

dv Vg - vi 

(^)v = ^ (A8) 
where the hnearization relation between entropy and volume of (B21) is used. Eq. (A5) 
yields 



(AT) 



dT 

This is just the Clapeyron equation. 



L = T—{vg-vi) (A9) 



APPENDIX B: COEXISTENT GAS-LIQUID SYSTEM 

In this part, we will debate thermodynamics and statistical mechanics of a coexistent 
gas-liquid system. 

Considering a coexistent gas-liquid system in equilibrium with particle number = 
Ng + Ni, volume V = Vg + Vi. For a given temperature T < Tc, densities of vapor and liquid, 
Ug and ni, are determined by Eqs. (16) and (17). First, we examine its thermodynamics. 
We have 

dUg = TdSg - PdVg + jjd N g (Bl) 

dUi = TdSi - PdVi + i^dNi (B2) 

So 

dU = TdS - PdV (B3) 
where U = Ug + Ui, S = Sg + Si , and dNg = —dNi. In the same way, we have 

dF = -SdT - PdV (B4) 

dG = -SdT + VdP (B5) 

Now we discuss statistical mechanics for VDW fluid. For the vapor and liquid parts, the 
partition functions are 

/ \ — TV 

1 (2TimkBTY 'aaN„n„nr l.^N, 



Z, = j^[ ) e/^-^"^ {Vg - Ngbf^ (B6) 



15 



respectively. We have 



(B7) 



Pg = ksT 



dVg I -Ugh ' 



fB8) 



(B9) 



(BIO) 



The partition function for the system is 



(Bll) 



Z — Zg ■ Zl 



(B12) 



We then define pressure and chemical potential for the coexisting system 



d 

P = knT [ —\nZ 



ngniksT 

,oV J J. ^ rii — Ug 



In ( — - 6 1 - In f — - 6 



n„ 



ni 



(B13) 



d 



knT 



knT 



T,V 

{ug + ni) + 



-ni \ng 

2'KmkBT 



n„ - ni \ni 



(B14) 



Infl-M+l 



It is easy to prove P = Pg = Pg and fi = fig = fii if the system is in equilibrium, i.e. when 
we have fig = and Pg = Pi. To do this, we only need to calculate 



2P-Pg- Pi 



Ug + ni 
m - n„ 



Pg -Pi) = 



(B15) 



2fl- fl -fll = g 

^ ni-Ug ^ 



{Pg -Pi) = o 



(B16) 



where = is used. 
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Under the conditions /x^ = /i; and Pg = Pi for a given temperature T. partition function 
Z = Zg ■ Zi has only one independent variable. Let x = Ng/N, then partition function can 
be written as 

^ ' ~ {xN)\[{l-x)N]\ \ ) ^ ' 

V% J \ J 

Next we will derive a thermodynamic relation of a coexistent gas-liquid system in equi- 
librium with fixed T and P. The molar entropy of the system is 

s = xsg + {l- x)si (B18) 

where Sg, si are molar entropy of gas and liquid respectively. From the relation 

V = xVg + (1 — x)vi (B19) 

we have 

a; = ^^,l-x = ^i^ (B20) 

Vg -Vl Vg- Vl 

substituting Eq. (B20) to (B18), then 

Sg- Sl VgSl - ViSg 

s = — V + — (B21) 

Vg - Vl Vg- Vl 

when T,P are fixed, Sg,si,Vg,vi are all constant. So from this relation, we can see that 
entropy is the linear function of volume under the given conditions. 
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